Shell-model calculations of stellar weak interaction rates: 
II. Weak rates for nuclei in the mass range A = 45 — 65 in supernovae environment 
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Based on large-scale shell model calculations we have determined the electron capture, positron 
capture and beta-decay rates on more than 100 nuclei in the mass range A — 45-65. The rates 
are given for densities pYe = 10^-10^" mol/cm'' and temperatures T — 10^-10^" K and hence are 
relevant for both types of supernovae (Type la and Type II). The shell model electron capture rates 
are significantly smaller than currently assumed. For proton-to-baryon ratios Ve = 0.42-0.46 mol/g, 
the beta-decay rates are faster than the electron capture rates during the core collapse of a massive 
star. 



PACS numbers: 26.50.+X, 23.40.-s, 21.60. Cs 
I. INTRODUCTION 

Astrophysical environments can reach very high densi- 
ties and temperatures. Under these conditions (temper- 
atures T larger than a few 10^ K), reactions mediated 
by the strong and electromagnetic force are in chemical 
equilibrium and the matter composition is given by nu- 
clear statistical equilibrium i.e. it is determined 
mainly due to the nuclear binding energies subject to the 
constraint that the total number of protons in the com- 
position balances the number of electrons present in the 
environment. Introducing the electron-to-baryon ratio 
(units of mol/g), this constraint can be formulated as 



(1) 



where the sum is over all nuclear species present and 
Zk, Ak, and Xk are the proton number, mass number 
and mass fraction of species k, respectively. Importantly, 
in these astrophysical environments the relevant density 
and time scales are often such that the neutrinos are ra- 
diated away, so that reactions mediated by the weak in- 
teraction are not in equilibrium. Thus, weak interaction 
rates play a decisive role in these environments changing 
Ye and hence the composition of the matter. 

Among these astrophysical environments are the two 
major contributors to the element production in the uni- 
verse: supernovae of type la and type II (e.g. |^). A 
type la supernova is usually associated with a thermonu- 
clear explosion on an accreting white dwarf. Hydrogen 
mass flow from the companion star in the binary system 
at rather high rates leads to steady hydrogen and helium 
burning on the surface, increasing the carbon and oxygen 
mass of the white dwarf. Finally carbon is ignited in the 
center of the star leading to a thermonuclear runaway. 
A burning front then moves outwards through the star 
at subsonic speed, finally leading to a detonation which 
explodes the star. Several issues (including the masses 
of the stars in the binary, the mass accretion history and 
composition, the matter transport during the explosion, 
the speed of the burning front) are currently still under 



debate (e.g. Q). It appears, however, established that 
electron capture will occur in the burning front driving 
the matter to larger neutron excess. As the observed iron 
in the universe is made roughly in equal amounts by type 
la and type II supernovae, type la supernovae should pro- 
duce the relative abundances within isotope chains in the 
iron mass region in agreement within a factor of 2 with 
the observed solar abundances |^ . Provided the electron 
capture rates are sufficiently well known (matter is in nu- 
clear statistical equilibrium inside a type la supernova), 
this requirement allows to put severe constraints on type 
la supernova models 

A type II supernova is related to the core collapse of 
a massive star. Here the core of a massive star becomes 
dynamically unstable when it exhausts its nuclear fuel. 
If the core mass exceeds the appropriate Chandrasekhar 
mass, electron degeneracy pressure cannot longer stabi- 
lize its center and it collapses. In the initial stage of the 
collapse, electrons are captured by nuclei in the nickel 
mass range, thus reducing Yf,. Associated is a decrease 
in degeneracy pressure and energy, as the neutrinos can 
still leave the star; both effects accelerate the collapse. 
With decreasing Y^, i.e. with increasing neutron excess 
of the nuclei present, P decay becomes more important 
and can compete with electron capture. 

Under the stellar conditions discussed above, the weak 
interaction rates are dominated by Gamow- Teller (GT) 
and, if applicable, by Fermi transitions. Bethe et al. M] 
recognized the importance of the collective GT resonance 
for stellar electron capture. Shortly after. Fuller, Fowler 
and Newman (usually abbreviated as FFN ||^-|ll|]) esti- 
mated the stellar electron capture and beta-decay rates 
systematically for nuclei in the mass range A = 45 — 60 
considering two distinct contributions. At first, these 
authors estimated the GT contributions to the rates 
by a parametrization based on the independent parti- 
cle model. The rate estimate has then been completed 
by including Fermi transitions and by experimental data 
for discrete transitions, whenever available. Unmeasured 
allowed GT transitions have been assigned an empirical 
value (log ft — 5). One of the important ideas in the sem- 
inal work by FFN was to recognize the role played by the 
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GT resonance in /5~ decay via the GT back resonance in 
the parent nucleus (the GT back resonance are the states 
reached by the strong GT transitions in the inverse pro- 
cess (electron capture) built on the ground and excited 
states, see allowing for a transition with a large 

nuclear matrix element and increased phase space. Until 
now the FFN weak interaction rates are key ingredients 
in supernova simulations. 

The general formalism to calculate weak interaction 
rates for stellar environment has been already given by 
Fuller et al. P-pl]|. What had not been possible at the 
time when FFN did their pioneering work was to solve 
the associated nuclear structure problem with the neces- 
sary accuracy and predictive power. Several years ago, 
Aufderheide, Mathews and collaborators |12|-[l4| pointed 
out that the interacting shell model is the method of 
choice for this job. In fact, using the newly developed 
shell model Monte Carlo techniques electron cap- 

ture rates for selected nuclei have been derived but 
it became apparent that shell model diagonalization cal- 
culations are preferable as they allow for detailed 
spectroscopy and do not have restrictions in their appli- 
cability to odd-odd and odd-A nuclei as the shell model 
Monte Carlo method has at low temperatures Be- 
fore calculating weak interaction rates for a large set of 
nuclei in the A = 45-65 mass region, it had to be proven 
that state-of-the-art shell model diagonalization is indeed 
capable of reliably solving the required nuclear struc- 
ture problems. This proof has been given in Ref. [ pO[ . 
This paper reported about large-scale shell model cal- 
culations covering the relevant mass range. The studies 
were performed at a truncation level in the pf shell at 
which the GT strength distributions are virtually con- 
verged. As residual interaction a slightly modified ver- 
sion of the wellknown KB3 interaction ||2l| has been used; 
the slight modifications correct for the small overbinding 
at the = 28 shell closure encountered with the origi- 
nal KB3 force |^^. In general, it has been demonstrated 
that the shell model reproduces all measured GT^ distri- 
butions very well and gives a very reasonable account of 
the experimentally known GT- distributions. Further, 
the lifetimes of the nuclei and the spectroscopy at low 
energies is simultaneously also described well. Ref. ||20| ] 
has therefore shown that modern shell model approaches 
have the necessary predictive power to reliably estimate 
stellar weak interaction rates. 

In this paper we will use the shell model approach 
of (20| and derive the weak interaction rates for more 
than 100 nuclei in the mass range A — 45-65 at a tem- 
perature and density regime relevant for supernova appli- 
cations. Our paper is organized as follows. In section 2 
we will repeat the derivation of the necessary formalism 
for the weak rates. The results are presented and ex- 
plored in section 3. In this section we will also compare 
them to the pioneering work of FFN. In this comparison 
we will find systematic differences. The origin for these 
differences will be explored and discussed in section 4. 



II. STELLAR WEAK INTERACTION RATES 
FORMALISM 

The definition of the stellar electron and positron cap- 
ture and /3-decay rates has been derived by Fuller, Fowler 
and Newman P~[ll[|. We will here repeat the formulae 
and ideas which will allow us to explain our strategy and 
procedure to calculate these rates and to point out dif- 
ferences with previous compilations. 



A. General formalism 

We computed rates for four processes mediated by the 
charged weak interaction: 

1. Electron capture (ec), 

{Z,A)+e- ^{Z-l,A) + iy. (2a) 



2. P+ decay (/3+), 

{Z,A) ^ {Z -l,A) + e+ 



(2b) 



3. Positron capture (pc), 

{Z,A)+e+ ^ {Z + 1,A) + D. (2c) 

4. /3" decay (/?"), 

{Z,A)^ {Z+l,A)+e- +D. (2d) 

The rate for these weak processes is given by 

ln2^ (2J. + l)e-W) 
^ -IT 4- G{Z,A,T) A.^^^^^.' (3) 

where the sums in i and j run over states in the parent 
and daugther nuclei respectively and the superscript a 
stands for ec, /?''", pc or The constant K is defined 
as 



K 



27r3(ln2)?i^ 



(4) 



where Gp is the Fermi coupling constant, Vud is the 
up-down element in the Cabibbo-Kobayashi-Maskawa 
quark-mixing matrix and = 1 is the weak vector 
coupling constant. K can be determined from superal- 
lowed Fermi transitions and we used K = 6146 ±6 s [ p2[ . 
G{Z,A,T) = Y.^^^v{-E^/{kT)) is the partition func- 
tion of the parent nucleus. Bij is the reduced transition 
probability of the nuclear transition. We will only con- 
sider Fermi and GT contributions which, however, has 
been shown to be quite sufficient: 



B,, ^ B,,{F) + B,,{GT). 



(5) 
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The GT matrix is given by: 



B,,{GT) 



9A 
9v 



011 E, 



.crHil 



2J, + l 



(6) 



where the matrix element is reduced with respect to the 
spin operator cr only (Racah convention and the 

sum runs over all nucleons. For the isospin rising and 
lowering operators, t± — [t^ ± iTy)/2, we use the con- 
vention = n; thus, '+' refers to electron capture 
and transitions and '— ' to positron capture and (3~ 
transitions. Finally, {gA/9v)ctf is the effective ratio of 
axial and vector coupling constants that takes into ac- 
count the observed quenching of the GT strength [p4[ . 
We use 



9v 



9A] ^Q.74(9± 



off 



gv 



(7) 



bare 



with (5yi/gy)baro = -1.2599(25) If the parent nu- 

cleus (with isospin T) has a neutron excess, then the 
GT_ operator can connect to states with isospin T — 1, 
T, T+1 in the daughter, while GT4. can only reach states 
with T+1. This isospin selection is one reason why the 
GT_|_ strength is more concentrated in the daughter nu- 
cleus (usually within a few MeV around the centroid of 
the GT resonance), while the GT_ is spread over 10- 
15 MeV in the daughter nucleus and is significantly more 
structured. 

The Fermi matrix element is given by: 



2J, + 1 



(8) 



In our calculations isospin is a good quantum number 
and the Fermi transition strength is concentrated in the 
isobaric analog state (IAS) of the parent state. Equa- 
tion M) reduces to. 



B,,{F)^T{T+1)-T,J,^, 



(9) 



where j denotes the IAS of the state i. We neglect the 
reduction in the overlap between nuclear wave functions 
due to isospin mixing which is estimated to be small (w 
0.5% H). 

The last factor in equation (||) , , is the phase space 
integral given by 



/3+ 



wp{Qij + wYF{Z, w) 

Se{w){\ - S^iQij + w))dw, 



wpiQ,j~w)^F{-Z + l,w) 



(10a) 



il-Sp{w))il-S,{Q,j-w))dw, (10b) 



wp{Q,j-wfF{Z + \,w) 

{l-Se{w)){l-S,{Q^,-w))dw, (10c) 



wpiQij + w)'^F{-Z, w) 

Sp{w){l - S,{Q^j + w))dw, (lOd) 



where w is the total, rest mass and kinetic, energy of the 
electron or positron in units of nieC^, and p = V^f^ — 1 is 
the momentum in units of mgC. We have introduced the 
total energy available in /3-decay, Qi, , in units of WeC^ 



Q^ 



1 



(11) 



where Mp , Md are the nuclear masses of the parent and 
daughter nucleus, respectively, while Ei,Ej are the ex- 
citation energies of the initial and final states. We have 
calculated the nuclear masses from the tabulated atomic 
masses neglecting atomic hinging energies, wi is the cap- 
ture threshold total energy, rest plus kinetic, in units 
of mgC^ for positron (or electron) capture. Depend- 
ing on the value of Qij in the corresponding electron 
(or positron) emission one has = 1 if Qij > — 1, 
or wi = \Qij\ if Qij < —1- SejSp, and S^, are the 
positron, electron, and neutrino (or antineutrino) distri- 
bution functions, respectively. For the stellar conditions 
we are interested in, electrons and positrons are well de- 
scribed by Fermi-Dirac distributions, with temperature 
T and chemical potential fi. For electrons. 



Se 



1 



exp 



fcT 



(12) 



with Ef, — 



The positron distribution is defined 



similarly with jip = —pLe. The chemical potential, ^le, 
determined from the density inverting the relation 



pYe 



1 



Sp)p^dp, (13) 



where Na is Avagadro's number. Note that the density 
of electron-positron pairs has been removed in (|l^ ) by 
forming the difference Se — Sp. 

In supernovae weak interactions with nuclei with mass 
numbers A — 45-65 occur at such densities that the neu- 
trinos can leave the star unhindered. Thus, there is no 
neutrino blocking of the phase space, i.e. Si, — 0. 

The remaining factor appearing in the phase space in- 
tegrals is the Fermi function, F{Z,w), that corrects the 
phase space integral for the Coulomb distortion of the 
electron or positron wave function near the nucleus. It 
can be approximated by 



F{Z,w)^2{l+j){2pR)-^<^'-r)^^!^±^e-y, (14) 



1 



where 7 = -^/l — (aZ)^, y — aZw/p, a is the fine struc- 
ture constant, and R is the nuclear radius. 

Finally, the calculation of the rates reduces to the eval- 
uation of the nuclear transition matrix elements for the 
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GT operator. The problem obviously lies in the fact that 
many states (can) contribute to the two sums over 
At first, the finite temperature allows the thermal popu- 
lation of excited states in the parent. Each of these states 
is then connected to many levels in the daughter nucleus 
by the GT operators. A state- by-state evaluation of both 
sums is still beyond present-day computer abilities. Be- 
fore we summarize our strategy to approximate the sums 
we recall that previous compilations of the stellar weak 
rates employed the so-called Brink hypothesis (e.g. [p8|): 
Let So{E) be the GT distribution in the daughter nu- 
cleus build on the ground state, then it is assumed that 
the distribution Si{E) build on the excited state in the 
parent at excitation energy Ei is the same as 5*0, but 
shifted in energy by E,, i.e. So{E) = S^{E + E,). This 
hypothesis has been tested in various shell model calcu- 
lations and is found to be valid for the gross structure of 
the GT distribution. However, it can be badly violated 
for specific transitions as they occur at low excitation 
energies (and are important for nuclear lifetimes). This 
tells us that Brink's hypothesis should not be employed 
if specific low-lying transitions (which usually exhaust a 
very small fraction of the total GT strength) dominate 
the rates or are important. Such situations may occur 
at low temperatures and densities. When the tempera- 
tures and densities are higher so that many states (a few 
tens or more) contribute, variations in low- lying transi- 
tion strengths tend to cancel and Brink's hypothesis be- 
comes a valid approximation. We will demonstrate this 
in detail later in this paper. 

B. Evaluation of the rates and checks 

Our strategy to calculate the weak rates (^) is best 
explained by considering a pair of nuclei {Z, A) and 
{Z + 1,A) connected by the weak processes under dis- 
cussion here. We have then calculated the GT+ distri- 
butions for all individual levels in the nucleus {Z + 1,A) 
at modest excitation energies and the GT_ strength dis- 
tributions for the low- lying levels in the nucleus {Z,A). 
Whenever experimental information about excitation en- 
ergies or GT transition strengths is available, the shell 
model results have been replaced by data. Otherwise the 
shell model predictions are used. The quality of these cal- 
culations is demonstrated in |Q where calculated spectra 
and lifetimes are compared to data. For even-even par- 
ents we considered explicitly the lowest O"*", 4+ and the 
two lowest 2+ levels which describe the spectrum typ- 
ically upto excitation energies of 2 MeV. Depending on 
the level spectrum, we adopted between 4 and 12 individ- 
ual states for odd-^ and odd-odd parent nuclei including 
at least the shell model spectrum at excitation energies 
below 1 MeV explicitly. As our shell model GT distri- 
butions have been calculated with 33 Lanczos iterations 
(see H^), a state in the parent nucleus is connected to 
33 ( for angular momenta J = 0) til 99 (for angular mo- 



menta J > 1) states in the daughter nucleus by the GT+ 
operator, while this magnifold is tripled for the GT_ op- 
erator due to the different isospin final states. (However, 
the T-l-1 component does not play a role in calculating 
the rates under the temperature/density conditions we 
are concerned with and is often omitted.) 

While the explicit consideration of the low-lying states 
guarantees a reliable description of the rates at low tem- 
peratures/densities, where individual transitions are of- 
ten decisive, states at higher excitation energies become 
increasingly important at the higher temperature and 
density regimes under consideration here. This is par- 
ticularly true for the P~ decay rates which is often 
dominated by the back-resonances under these condi- 
tions We have therefore supplemented the con- 
tribution of the low-lyin states by the back-resonance 
contributions which can be derived from the low-lying 
contributions in the inverse process, i.e. the states pop- 
ulated by electron capture on {/3~ decay of) low-lying 
states became the back-resonances in the /3~ -decay (elec- 
tron capture). The back-resonance contribution defined 
this way does not exhaust the total GT strength built 
on this excited state. In particular, the capture of high 
energy electrons in the back-resonance states will lead 
to states in the daughter which are not included in our 
model. To correct for these missing transitions we have 
employed the Brink hypothesis, i.e. we calculate the total 
GT strength and the centroid Ec for the parent ground 
state and place the strength of the back-resonance state 
at an energy Ei + Ec, where Ei is the energy of the back- 
resonance state in the parent. 

Fig. |l| illustrates the various contributions to the two 
sums in the rate formulae. In evaluating these sums we 
typically consider several hundred states in both the par- 
ent and the daughter nucleus. The partition function is 
consistently derived from the same parent states. We do 
not introduce a cut-off of levels at the particle separation 
thresholds [||J2|]. 

To demonstrate the convergence of our rates we will ex- 
plicitly discuss the pair of nuclei ®^Co and ^^Ni. Due to 
Rcf. [g^, ^-^Co is among the most important /3-decaying 
nuclei in a presupernova collapse at densities around 10® 
g/cm'^, while electron capture on ^'^Ni plays a moder- 
ate role at the same astrophysical conditions. The Qp 
value for this pair of nuclei is 3.672(20) MeV. Thus un- 
der laboratory conditions ^'^Co /3-decays to ^■^Ni, leading 
dominantly to the excited 5/2 state at 87 keV with an 
log ft value of 4.8(1). Our shell model calculation agrees 
with the data. 

Our calculation of the stellar electron capture (and /3+ 
decay) rate explicitly considers the two lowest 1/2^ and 
3/2~ states and the lowest 5/2^ state in ^^^Ni, which 
comprises the experimentally known spectrum upto an 
excitation energy of 1 MeV. Furthermore we include the 
back-resonances from the inverse reaction which intro- 
duces a total of 6 states in the excitation energy interval 
between 1 and 2 MeV, while experimentally 8 states are 
known in this energy regime. Between 2 and 5 MeV, 



4 




(Z+1,A) 

FIG. 1. Sketch of the GT transitions considered in the 
present estimate of the stellar weak interaction rates. Our 
evaluation of the sums in Eq. ^ explicitly includes the GT_ 
and GT-i- strength distributions for the few lowest states in a 
nucleus. This is supplemented by the back-resonances which 
are derived from the GT strength distributions of the in- 
verse process, as described in the text. Experimental data 
for the energies and GT transition strengths have been used, 
whenever available. Otherwise the shell model results are 
adopted. Finally the Brink hypothesis (see text) is employed 
to derive the remaining GT strength and centroid position 
for back-resonance states. In the sketch these transitions are 
shown by dashed lines, while the GT transitions, for which 
either data or shell model results are used, are drawn by solid 
lines. 



we include 16 more states which, however, are probably 
not fully converged within our Lanczos procedure and 
thus represent 'averaged GT states' rather than physical 
states. 

Table | demonstrates the convergence of the stellar 
^^Ni electron capture rate with the number of initial 
states. We have calculated the rates for various densities 
between pYe = 10^ and 10^° mol/cm^. The temperatures 
have been chosen from the FFN temperature/density 
grid to be the closest to the expected stellar trajectory 
in a supernova collapse |Q . The calculations have been 
performed considering only the ground state, the lowest 
2 and 3 states, all 5 states for which shell model GT+ dis- 
tributions have been explicitly calculated and finally for 
our full procedure including individual states and back- 
resonances. 

The calculation performed by restricting the sum over 
initial states to only the ground state resembles the appli- 
cation of 'Brink's shift hypothesis', as with this assump- 
tion the nuclear matrix elements and the phase space fac- 
tors loose their dependence on the parent state and the 
sum over initial states cancels the partition function [ p8[ . 
(We point out that this calculation is, however, not the 
same as those in Refs. [Hl0[|2|] as we adopt the shell 



model GT_|_ strength distribution, while the other au- 
thors had used an empirical parametrization of the GT 
strength (see below).) As expected, the apphcation of 
the Brink hypothesis is not a good approximation at low 
temperatures and densities where the rate is sensitive 
to low-lying transitions which can vary strongly between 
various initial states. The assumption becomes, however, 
quite acceptable at higher densities where enough high- 
energy electrons are available to effectively reach the cen- 
troid of the GT distribution in the daughter. As stated 
above. Brink's shift hypothesis works well if centroids 
are compared. To demonstrate this quantitatively we 
have calculated the GT centroids {Egt) for the 5 initial 
states. A measure for the validity of the Brink hypothesis 
is then given by the difference Eqt — Ex , where Ex is the 
shell model excitation energy of the initial state. If mea- 
sured relative to the ^^Ni ground state, this difference is 
13.7 MeV for the ground state, while we find 13.7 MeV 
(5/2-), 13.7 MeV (3/2"), 13.9 MeV (excited 3/2") and 
14.1 MeV (excited 1/2") for the other states. 

In Table I we also observe that, even at low densities, a 
convergence is rather fast achieved if the rate comprises 
a thermal average over a few states. 

The ^^Co /3-decay rate has been calculated by ex- 
plicitly considering the GT_ distributions for the low- 
est 7/2" (ground state), 3/2" (at 0.995 MeV), 5/2" (at 
1.437 MeV) and 1/2" (at 1.889 MeV) states. Only the 
1/2^ state has a definitive spin assignment, the spin for 
the other states has been assigned based on our calcula- 
tions. These contributions are then supplemented by the 
back-resonances obtained from the GT+ distributions of 
the inverse reaction. These include another 4 states in 
the excitation energy interval between 1 and 2 MeV and 
18 (average) states between 2 and 5 MeV. Due to their 
construction, the back-resonances have only transitions 
to the 5 lowest levels in ^'^Ni. The missing decay chan- 
nels of these states are, however, significantly suppressed 
due to the strong Qfj dependence. 

In Table || we study the convergence of the ^^Co /3~ 
decay rate. The full rate is compared to rates in which 
the contributions from the 4 individual levels are kept, 
but the back-resonance contributions are cut at 5 MeV 
and 3 MeV in the parent nucleus. Additionally we have 
calculated the rates only from the GT_ distributions of 
the lowest 4, 2, and 1 states, totally neglecting the back- 
resonances. 

We find that at low temperature/densities the rate can 
be calculated solely from the individual levels in a good 
approximation; even the decay rate determined from the 
ground state alone is already a fair approximation. At 
moderate and higher densities the back-resonances be- 
come increasingly more important. This is easily under- 
stood by the fact that at these high densities the Fermi 
energy of the electrons gets so high that decays from low- 
lying states are effectively blocked. Correspondingly the 
(3 decay rates decrease with increasing density. But we 
also note from Table || that the rate still converges rather 
rapidly and that considering the back-resonances upto 5 
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TABLE I. Convergence of the stellar Ni electron capture rate, logA"^*^, along the stellar trajectory during a supernova 
collapse (defined by temperature Tg and density p). Our full rate is compared to rates in which only the lowest 5, 3, 2, and 
1 states in ®^Ni are considered. The 'no Brink' calculation has been performed neglecting those rate contributions which are 
shown by dashed arrow and are labelled 'Brink hypothesis' in Fig. nl. 



Tg logpVj full no Brink 5 states 3 states 2 states 1 state 

1 7 -8.46 -8.46 -8.62 -8.59 -8.46 -9.62 

4 8 -5.18 -5.18 -5.29 -5.25 -5.12 -5.83 

5 9 -1.74 -1.77 -1.77 -1.72 -1.61 -1.93 
7 10 1.78 1.63 1.74 1.78 1.80 1.84 



TABLE IL Convergence of the stellar ^^Co /3~ decay rate, log A'' , along the stellar trajectory during a supernova collapse 
(defined by temperature Tg and density pr). Our full rate is compared to rates in which back-resonant states are only considered 
upto 5 MeV and 3 MeV, respectively, and to rates which are calculated from the GT- distributions of the lowest 4, 2, and 1 
states in "^Co. 



Tg 


log pYe 


full 


cut 5 MeV 


cut 3 MeV 


4 states 


2 states 


1 state 


3 


7 


-1.60 


-1.60 


-1.60 


-1.61 


-1.62 


-1.64 


4 


8 


-1.77 


-1.77 


-1.79 


-1.89 


-1.90 


-1.97 


5 


9 


-2.55 


-2.55 


-2.71 


-3.52 


-3.65 


-4.03 


7 


10 


-5.43 


-5.51 


-5.99 


-7.12 


-7.30 


-7.71 



MeV excitation energy gives a sufficient approximation. 

We mention that, due to at high densities nuclei 
in the supernova environment are expected to be more 
neutron-rich than ^■^Co. Thus the relevant Qfj values 
for (3 decay increase making the final-state blocking by 
electrons less effective and reducing the importance of 
back-resonant states relative to the low-lying individual 
states which we have considered. 



III. STELLAR WEAK RATES 

We have calculated the stellar weak interaction rates 
(electron and positron capture, (3~ and /?"'" decay) for 
more than 100 nuclei in the mass range A = 45-65. The 
rates have been calculated for the same temperature and 
density grid as the standard FFN compilations . 
An electronic table of our rates is available from the 
present authors upon request. We have also prepared 
a table in which the electron capture rates are presented 
in terms of the 'effective rates' introduced and defined by 
Fuller et al. [ pl| which allow a more reliable interpolation. 

Examples of our rates are shown in figures ^ and ^ 
using the pairs of nuclei (^^Co, ^^Ni) and (^^Fe,^^Co) 
as typical examples. The latter pair includes an even- 
even and odd-odd nucleus. The Qp value for these nuclei 
is negative and in the laboratory ^^Co decays to ^^Fe 
by electron capture. To calculate the electron capture 
rate we have considered the 5 states in ^^Co below 1 
MeV excitation energy and the lowest 1+ state (at 1.72 
MeV) explicitly, supplemented by the back-resonances. 



The ^^Fe (3 decay rate can be calculated from the back- 
resonances due to the negative Q value. 

The electron capture rates on both nuclei ^"^Ni and 
^^Co increase with temperature and density. Due to the 
negative value, electron capture is possible from all 
states in ^^Co and thus the increase with temperature is 
rather mild. However, in the case of ^'^Ni electron cap- 
ture is quite sensitive to temperature (at low densities), 
as it has to overcome a threshold of nearly 4 MeV which 
at the low densities (pYg — 10^ and 10^° g/cm'^) is mainly 
achieved by thermal population of excited states in the 
parent nucleus, as the Fermi energies of the electrons 
(1.2 MeV and 2.4 MeV, respectively for Tg = 1) is still 
noticeably smaller than the Qp value of 3.67 MeV. How- 
ever, at the higher densities electrons with energies above 
Ql3 are sufficiently available to allow for efficient capture. 
In this case, the rate becomes only slightly dependent on 
temperature. 

The '^^Fe /3-decay rate shows a very steep tempera- 
ture dependence and decreases with density. Both ef- 
fects are readily explained. Due to the threshold of about 
4.5 MeV, (3 decay is only possible from moderately ex- 
cited states which have to be populated thermally lead- 
ing to the strong temperature dependence. As the elec- 
tron Fermi energy increases with density, an increasingly 
larger part of the phase space gets Pauli-blocked. Conse- 
quently the /3-decay rate decreases with increasing den- 
sity. Note that the centroid of the GT_|_ strength distribu- 
tions of the ^^Co ground state is at an excitation energy 
of around 8 MeV in ^^Fe. Thus at densities larger than 
pYg — 10° mol/cm'^ (corresponding to an electron Fermi 
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energy of about 5 MeV for Tg = 5), the strong transitions 
associated with the centroid of the back-resonances are 
getting Pauh-blocked explaining the larger decrease with 
density even at high temperatures. 

^•^Co can (3 decay from all states and indeed already 
the ground state has a strong GT transition to the ex- 
cited state in ^^Ni at 87 keV. Furthermore, the first ex- 
cited state is at nearly 1 MeV excitation energy. Conse- 
quently, the ^'^Co /3-decay rate shows only a mild tem- 
perature dependence at low densities. Pauli-blocking by 
the electrons in the final state becomes, however, impor- 
tant at higher densities introducing effectively a thresh- 
old for the /3 decay which has to be overcome by thermal 
population of excited states in the parent. As a result, 
the rate develops an increasing sensitivity to tempera- 
ture. The centroid of the GT+ strength distribution of 
the ^^Ni ground state is at an excitation energy of around 
2.5 MeV in ^^Co. Consequently Pauli blocking by elec- 
trons becomes increasingly important for densities above 
10^ g/cm^. 

The positron distribution does not play a role for the 
f3~^ decay rate, which is virtually independent on density 
for both nuclei ^^Ni and ^^Co (this is already apparent 
from the FFN rates [|,|o)). The ^^Co /3+ rate depends 
only mildly on temperature. This is different for ^^Ni 
where a threshold has to be overcome by thermal popu- 
lation in the parent. For both nuclei the f3~^ decay rate is 
noticeably smaller than the electron capture rate which 
generally dominates the weak rates for charge-decreasing 



nuclear transitions under supernova conditions. 

The positron capture rate decreases with density, but 
increases with temperature. Both dependencies are 
caused by the positron distribution where an increasing 
number of high-energy positrons gets available by rais- 
ing the temperature or lowering the density. The lat- 
ter is caused by the fact that the degeneracy parameter 
fi/kT is negative for positrons. Due to its threshold for 
positron capture, the rate on ^^Fe shows the steeper tem- 
perature dependence. Again, the positron capture rates 
are usually smaller than the competing f3~ rates under 
supernova conditions for the nuclei of interest here. 

The most interesting question clearly is: How do the 
shell model rates compare to the FFN rates? 

To answer this question, we have calculated the ratio 
Affn/Asm, where Affn and Asm are the FFN and shell 
model rates, respectively. The FFN rates are taken from 
the electronic file available at |30|. For the comparison 
we choose 4 different temperature and density grid points 
(r9,log(pre)): (3,7), (5,8), (5,9) and (10,10). Ratios for 
the two important weak processes, electron capture and 
(3~ decay, are plotted in figures |^ and || for 4 different 
chains of isotones equally spanning the mass range be- 
tween A = 50-60. 

To understand the differences observed in figures ^ 
and H we have to recall how the FFN rates have been 
derived. At first, these authors considered experimental 
data for discrete transitions, whenever available (like we 
do in the present shell model rates). The main contribu- 
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FIG. 4. Ratio of the FFN and shell model electron capture 
rates for nuclei in the mass chains A = 51, 54, 57, and 60 
for selected values of temperature and density, T and density, 
log{pYe). The charge numbers refer to the parent nucleus. 
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FIG. 5. Ratio of the FFN and shell model /3-decay rates for 
nuclei in the mass chains A — 51, 54, 57, and 60 for selected 
values of temperature T and density, log(pYe). The charge 
numbers refer to the daughter nuclei. 



tion to the FFN rates usually comes from the so-called 
GT resonance which they parametrized on the basis of 
the independent particle model and which represents the 
total GT strength by a single state. Here it is not so im- 
portant that the authors did not explicitly consider the 
quenching of the GT strength with respect to the inde- 
pendent particle modelQ. More relevant is where FFN 
placed the GT resonance. Here we will concentrate on 
the electron capture. The differences in the (3~ decay 
rates follow then from a discussion of the back-resonance 
contributions. 

FFN estimated the GT resonance energy Ectr from 
3 distinct contributions: 

EgTR — ^Esp + lS.Eph + AEpair- (15) 

The single particle term is calculated as follows. Start- 
ing with the independent particle model wave functions 
for the ground states of the parent and daughter nucleus, 
protons are acted on with the GT+ operator leading to 
final neutron states, which, within the independent parti- 



*The quenching in the GT+ strength had not been estab- 
lished at the time FFN calculated the rates. In later discus- 
sions [|ll| these authors point to this effect and showed a way 
how to incorporate quenching effectively into the rates. 



cle model, correspond to Ip-lh excitations of the ground 
state (or is the ground state). The corresponding exci- 
tation energy is readily calculated, where FFN used the 
single particle energies of Seeger Q . If the GT+ opera- 
tion can lead to several different final states the excitation 
energy of the resonance is computed taking the weighted 
average of the different transitions. AEph is the particle- 
hole repulsion energy which has to be supplied to pull a 
neutron out of the daughter ground state. For simplicity, 
FFN put AEph = 2 MeV for aU nuclei. Finally, FFN 
argued that there is a penalty energy which has to be 
paid to break a neutron pair if there is an even number 
of neutrons in the daughter ground state. When appli- 
cable, this pairing energy is approximated by 2 MeV for 
all nuclei. 

In previous publications we have already pointed 
out that the shell model rates for important electron- 
capturing nuclei along the stellar trajectory (due to the 
ranking given in these are generally odd-A and odd- 
odd nuclei) are usually smaller than the FFN rates. As 
one potential reason for this difference we have noted 
that the GT centroids for odd-A and odd-odd nuclei are 
usually at higher energies than assumed by FFN. We will 
discuss this difference and its potential origin in details 
below, but we note here that there are two other im- 
portant ingredients (low lying transitions and Qfs values) 
which can lead to differences between the shell model 
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and FFN electron capture rates. In fact, for ^-'^Ti and 
^^V we find from fig. ^ that the shell model rates are 
larger than the FFN rates at low temperature/density. 
For these cases the difference is due to the fact that the 
shell model predicts low-lying strength for the ground 
states (in ^^V to the second excited 5/2 state in ^^Ti, 
in ^^Ti to the lowest 1/2 state in ^^Sc) which are larger 
than the standard assignment (log ft = 5) used in |^-[ic|] 
for experimentally not known transitions. At larger tem- 
peratures/densities the low- lying strengths becomes less 
important as the capture proceeds mainly to the GT res- 
onance; consequently the FFN rates are then larger than 
the shell model rates due to differences in the position of 
the centroid. (Note that the GT+ strength distribution 
for ^^V has been measured and it agrees nicely with the 
shell model results pof.) 

The nuclei ^^Mn and ^°Fe serve as examples for an- 
other source of differences between the shell model and 
FFN rates; the latter being noticeably smaller than the 
shell model rates at low temperatures/densities. This 
is due to the use of different values. FFN had to 
rely on the systematics available at the time, while mod- 
ern compilations indicate that the values for these 
nuclei are about 950 keV (^'^Mn) and 1.7 MeV (^°Fe) 
smaller than adopted by FFN. Obviously the too large 
Q value suppressed the electron capture at low temper- 
atures/densities. With increasing temperature/density 
the electron Fermi energy grows strongly reducing the 
sensitivity to differences in the Q value. 

More generally, one expects that the electron cap- 
ture rates become less dependent on details of the GT 
strength distributions with increasing electron Fermi en- 
ergies. This explains why the deviations between the 
FFN and shell model rates reduce with increasing tem- 
perature/density. In this limit, the FFN rates should 
be still slightly larger than the shell model rates due to 
the neglect of the quenching of the total GT strength in 

From the above discussion about the various contri- 
butions to the rates, one might distinguish 3 different 
temperature/density regimes along a stellar trajectory. 
At low (T, p) specific low-lying transitions can be quite 
important, supplementing the rate contribution from the 
GT resonance. This is particularly the case for electron 
capture, if the Q value only allows capture of high-energy 
electrons from the tail of the Fermi-Dirac distribution. 
At intermediate (T, p) the rates are usually dominated 
by the strong transitions involving the GT resonance. At 
high (T, p) (when is large compared to Qij for transi- 
tions to the GT centroid) the rate becomes insensitive to 
the energy dependence of the GT distribution and hence 
the rate depends only on the total GT strength. 

Differences between the shell model results and the 
FFN assumptions in the 3 ingredients (GT centroid en- 
ergy, low-lying strength, Q values) lead also to differences 
in the decay rates, as is shown in figure |^. Although 
the FFN /3 decay rates are usually somewhat larger than 
the shell model rates, no general picture emerges in this 



comparison, as the rates are usually given by the sum 
of low-lying transitions and contributions from the GT_ 
distribution as well as of the back-resonances. As is dis- 
cussed in the misplacement of the GT^ centroid ef- 
fects the contribution of the back-resonances to the /3- 
decay rates. In particular, the back-resonances in odd- 
odd parents (the GT+ centroids of even-even nuclei have 
been often placed at too high energies by FFN) can 
be thermally excited more easily than assumed by FFN 
resulting in slightly larger /3-decay rates for these nu- 
clei. Examples are the odd-odd nuclei ^"Mn and ^°Cu. 
As FFN did not consider the quenching of the GT_|_ 
strength, the contribution of the back-resonances is some- 
what overestimated in FFN explaining the slightly larger 
FFN /3-decay rates for ^''V or ^''Mn. Differences in the 
adopted Q-values effect mainly neutron-rich nuclei. An 
example here is the /3-decay of ^^Cr, where FFN used the 
Q-value of 6.56 MeV rather than 5.60 MeV. 

The p/-shell nuclei discussed here have not too extreme 
Ye values and are therefore important at the earlier stages 
of the presupernova collapse involving low and intermedi- 
ate (T, p) values. As in these regimes the energy position 
of the GT centroid plays an essential role, we will now 
investigate in details the differences between the place- 
ment of this centroid in the FFN compilation with the 
shell model rates (and the data). 

In 1 32 1 we have pointed out that there are systematic 
differences between the placement of the GT centroid in 
the FFN compilations and the shell model results (and 
data, if available). It has been noted that the differences 
apparently depend on the pairing structure of the par- 
ent nucleus. For even-even parents, the GT+ centroid 
is calculated at slightly smaller energies than assumed 
by FFN, while it has been noticed that for odd-A and 
odd-odd nuclei the centroid is at higher excitation en- 
ergies than parametrized by FFN. The consequences for 
the electron capture and /3-decay rates are obvious. If the 
kinematics is such that electron capture is dominated by 
transitions to the GT resonance, the FFN rates should 
be larger than the shell model rates as, for given temper- 
ature and density, less electrons are available to capture 
to the centroid if it resides at higher energies. This effect 
is strongest for odd-odd nuclei where the differences be- 
tween the FFN and shell model rates can be larger than 
2 orders of magnitude, i.e. for ^°Co (fig. ^ which is usu- 
ally considered to be among the most effective electron- 
capturing nuclei along the collapse trajectory of a super- 
nova. 

Also for the /3-decay rates the ratios systematically de- 
pend on the pairing structure of the parent. Here, how- 
ever, the shell model rates are similar to the FFN rates 
for odd-odd parents (they are even often slightly larger), 
while they are smaller for odd-A and even-even nuclei. 
The reduction is usually largest for even-even nuclei. 

The systematics assumed by FFN for the GT reso- 
nance energy in the daughter is most easily illustrated 
for parent nuclei with only fj/2 protons and more than 
28 neutrons (so that the fj/2 orbital is blocked for GT+ 
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transitions). Then the single particle contribution to the 
GT resonance energy is about 1.8 MeV, if the neutron 
number is smaller A'' < 32, or 0, if > 32. Thus one 
finds EcTR ~ 3.8 (2.0) MeV for even-even parent nu- 
clei and 5.8 (4.0) MeV for odd-odd parents. For odd-A 
parents one has Eqtr = 3.8 (2.0) MeV if the neutron 
number in the parent is even and 5.8 (4.0) MeV, if it is 
odd. Here the numbers in parentheses refer to parent 
nuclei with A'' > 32 for which the P3/2 neutron orbitals 
are completely occupied. 
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FIG. 6. Comparison of the FFN GT+ resonance energies 
with the shell model centroid energies for various nuclei with 
Z < 28 and iV > 28. 



Figure || compares the shell model GT centroids with 
the FFN estimates of the GT resonance energy for about 
45 nuclei with proton numbers Z < 28 for which the 
above given systematic of the FFN estimates apply. In- 
deed one clearly observes that the FFN GT resonance 
energies cluster around 2, 4, and 6 MeV. Glearly one ex- 
pects that the shell model GT centroid energies are more 
scattered as the residual interaction fragments the GT 
strength and structure effects will also effect the GT dis- 
tributions. But a more striking result is found if we com- 
pare the GT centroids for even-even, odd-A and odd-odd 
parent nuclei separately. Here we find that, compared 
to the shell model centroids, FFN places the GT reso- 



nance energy usually at too high an energy for even-even 
parents and often at too low an energy for odd-odd par- 
ents. The situation is more interesting and also more 
telling for odd-A parents. Note that all those nuclei, 
for which FFN estimated the GT resonance energy by 
about 6 MeV, have parent nuclei with an odd number of 
neutrons and thus require an additional pairing energy 
in the FFN estimate. Gompared to the shell model GT 
centroids the FFN estimates for these nuclei are too high. 
Then compare the nuclei with an even neutron number 
in the parent. They do not acquire the additional pair- 
ing energy in the FFN estimate. Compared to the shell 
model estimate, the FFN GT resonance energies are usu- 
ally too low. 

TABLE III. Comparison of the experimental centroids of 
the GT resonance with the FFN GT resonance energies and 
the shell model centroids. The data are from [^-|3^. All 
energies are in MeV. 



nucleus 


data 


FFN 


shell model 


5^Fe 


3.7 ±0.2 


3.80 


3.78 


s^Fe 


2.6 ±0.2 


3.78 


2.60 




3.6 ±0.2 


3.76 


3.75 


60Ni 


2.4 ±0.3 


2.00 


2.88 


62Ni 


1.3 ±0.3 


2.00 


1.78 




0.8 ±0.3 


2.00 


0.50 


51y 


4.8 ±0.2 


3.83 


5.18 


^^Mn 


4.1 ±0.3 


3.79 


4.57 


58 Co 


4.4 ±0.3 


2.00 


5.05 



Thus it is obvious from this comparison that there is a 
different dependence on the pairing structure of the par- 
ent ground state between the FFN assumptions and the 
shell model results. But, which is correct? Experimen- 
tally the GT+ strength distribution has been studied for 
several even-even (^^-sepe, 58,60,62,64^^^^ 3 nu- 
clei (^^V, ^^Mn, ^^Go) in this mass range. As pointed 
out in the GT centroid for the even-even parents 
are generally at lower excitation energies in the daugh- 
ter than for the odd-^ nuclei. If the experimental cen- 
troids are compared to the FFN estimates one finds the 
same trend as discussed above for the comparison with 
the shell model centroids. Here we have calculated the 
experimental centroids from the measured GT distribu- 
tions upto 8 MeV. For ^^'^^Ni, however, we only consider 
the peak in the GT distribution corresponding to the 
GT resonance. The tail of the experimental distribution 
(see Fig. 1 in Q), if real, is most hkely due to states 
outside of our present model space. As demonstrated in 
Table [H, the FFN GT resonance energy is usually at 
too high an energy for even-even parents, while it is at 
too low an energy for the 3 odd-A nuclei, investigated 
experimentally. We note, however, that these 3 nuclei all 
have an even neutron number in the parent and just do 
not allow to explore the assumptions which FFN made 
concerning the pairing energy contribution to the GT res- 
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onance energy. We stress that our shell model centroids 
are generally in good agreement with the data, as in fact 
the shell model calculations reproduce all measured GT+ 
strength distributions quite satisfactory. 

The differences in placement of the GT centroids ex- 
plain the differences between the FFN and shell model 
weak interaction rates. As FFN assumes the GT reso- 
nance energy for odd-odd and odd-A parent nuclei with 
an even neutron number at lower energies than placed 
by the shell model, electron capture to these strong tran- 
sitions is easier and hence the FFN electron capture 
rate is significantly larger than the shell model rate for 
these nuclei. For electron capture on even-even nuclei 
and for odd-A nuclei with an odd neutron number, FFN 
have generally placed the GT resonance energy at lower 
daughter energies than predicted by the shell model (or 
the data). This effect alone would make the FFN rates 
smaller than the shell model rates. However, it is largely 
compensated by the fact that FFN did not consider the 
quenching of the GT strength and by their consideration 
of experimentally known transitions at low excitation en- 
ergies. Taken together, the FFN rates are also for these 
parent nuclei usually somewhat smaller than the shell 
model rates, with the notable exception of ^^Ni [|l8|. For 
reasons which will become apparent in the next section, 
the FFN rates on very neutron-rich odd-A nuclei with 
odd neutron number are also often smaller than the shell 
model rates. 

Remembering the importance of the back-resonances 
for the P~ decay rates, the effect of the different place- 
ments of the GT centroid is obvious. In even-even nuclei 
and odd-A nuclei with even neutron number the shell 
model studies place the back-resonances at higher exci- 
tation energies than assumed by FFN. Correspondingly, 
its thermal population becomes less likely and hence 
the contribution of the back-resonances to the de- 
cay rates decreases. On the contrary, experimental data 
and the shell model calculations indicate that the back- 
resonances reside actually at lower excitation energies in 
odd-odd nuclei than assumed by FFN. Consequently, the 
contribution of the back-resonances to the f3~ decay rate 
of odd-odd parent nuclei should be larger than assumed 
in the FFN rates, which is indeed the fact for several 
nuclei like '^'^'^^Mn and ^^Co. The effect of the misplace- 
ment of the GT centroids on the /3~ rates has already 
been discussed in |l^,|lj,^ . 

In the next section we will argue why the systemat- 
ics of the shell model GT centroids is correct. Further- 
more, in connection with the well-established systemat- 
ics of the GT centroids for the GT_ transitions, we will 
give a rather simple parametrization for these quantities. 
Although these arguments still have to be delivered, we 
will close this section pointing out that the present shell 
model rates are more reliable than the FFN ones and in 
fact should represent a fairly accurate description of the 
nuclear structure problem required to derive these rates. 
Then the weak interaction rates adopted in supernova 
simulations should be revised. We will speculate briefly 



below about possible consequences which the shell model 
rates might have for the presupernova core collapse. 



IV. SYSTEMATICS OF THE GT CENTROID 
ENERGIES 

We begin our discussion by recalling that the system- 
atics of the GT centroid energy Ear is well understood 
for GT_ transitions. In fact, one finds in a very good 
approximation |33] 



Egt - El AS = a + b- 



(N-Z) 
A ' 



(16) 



The constants a, b can be derived by fit to measured GT_ 
strengths p9]. 




FIG. 7. The shell model GT_ centroid energy, relative to 
the energy of the isobaric analog state (upper panel) and to 
the daugther ground state (lower panel), for the Mn, Fe, Co, 
and Ni isotopes as function of (A'^ — Z)/A. The pairing en- 
ergy in the mass splitting between parent and daughter nuclei 
shifts the GT_ centroid energy for odd-A and odd-odd par- 
ents up in energy with respect to even-even parent nuclei. 



In Fig. ^ we demonstrate that the centroids of our 
shell model GT_ strength distributions indeed exhibit 
this {N — Z) I A dependence. Our results show a spread 
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of about 1 MeV which can be interpreted as the uncer- 
tainty introduced by nuclear structure effects. In passing 
we note that for nuclei with large neutron excess the GT_ 
centroid energy can be below the IAS energy. As the fig- 
ure is compiled for the isotope chains with Z = 25 — 28 it 
comprises all 3 different types of parent nuclei: even-even, 
odd-odd and odd-A. Clearly no systematic dependence 
on the pairing structure is found. Obviously such a de- 
pendence shows up if the GT centroid energies are plot- 
ted as measured relatively to the daughter ground state 
energy (lower panel of figure |^) . It is trivially introduced 
by the differences in pairing energy between the various 
parent and daughter nuclei. As a consequence GT_ cen- 
troids form now 3 distinguished bands: one for even-even 
nuclei, one for odd-^ nuclei, and one for odd-odd nuclei 
the latter two shifted up in energy by about twice or four 
times the pairing energy, respectively. 

As we will show now the same behavior is found for 
the centroids of the GT+ strength distributions. Before 
we do so, however, it is useful to recall the origin of the 
(A'' — Z)/A dependence of the GT_ centroids. Here we 
will follow the nice discussion given by Bertsch and Es- 
bensen |40|. Then the average excitation energy of the 
GT operator is given by 



GT 



dEES{E) 



(a-t^ [H, (Tt±]) 



(17) 



where the expression is valid for both GT_ and GT+ op- 
erators. Several pieces of the Hamiltonian do not com- 
mute with the GT operator. While the most of these 
contributions are cancelled in building the difference with 
the IAS energy, the Wo-rf i • o'2Ti ■ T2 residual interaction 
gives rise to an energy shift that in the Tamm-Dancoff 
approximation is given by ||40| 



AE^r 



3A 



(18) 



where v^rP is the product of the integrated strength and 
the averaged ground state density. (In building the dif- 
ference with the IAS energy a similar contribution arising 
from the VrTi •T2 interaction has to be subtracted.) Sp is 
the total GT strength, which for the GT_ operator and 
neutron-rich nuclei, can be approximated by the Ikeda 
sumrule, 5*^- = 3(A^ — Z). Upon substitution into (18), 
one finds the desired {N — Z) /A dependence. 

For the GT+ operator one derives at the same ex- 
pression for lS.E„T, however, now considering Sp+. If 
one measures the GT+ centroid from the parent ground 
state, the contribution from the other terms of the Hamil- 
tonian are largely cancelled and one has upto a con- 
stant (reflecting for example the spin-orbit splitting) 
Egt = const -|- ErjT- The problem just reduces to find an 
appropriate parametrization of the total GT+ strength. 
As has been pointed out in Ref. , the presently avail- 
able data for pf shell nuclei suggest a scaling of the total 
GT+ strength like 



(20- A^,) 



(19) 



where Z^, are the number of valence protons and neu- 
trons in the pf shell, respectively. This dependence cor- 
responds to a generalized BCS model with pure proton 
and neutron pairing |41| . 
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FIG. 



8. Upper panel: shell model GT+ centroid energy, 
relative to the parent ground state energy, for various nuclei 
in the mass range A = 50 — 60 as function of Z„(20 — Nv)/A. 
Lower panel: shell model GT+ centroid energy, relative to the 
daughter ground state energy, for the same nuclei as function 
of (TV — Z)/A. The pairing energy in the mass splitting be- 
tween parent and daughter nuclei shifts the GT+ centroid 
energy for odd-y4 and odd-odd parents up in energy with re- 
spect to even-even parent nuclei. 



We have calculated the GT+ centroids for about 75 
nuclei in the pf shell and if we plot these energies with 
respect to the parent ground state energies, they closely 
follow the Zy ■ (20 — A'^^„) /A rule. Importantly, like for the 
case of the GT_ centroids, no dependence on the pairing 
structure of the parent nucleus is observed. However, 
such a dependence, as we have seen above for the GT_ 
centroids, is introduced if one measures the GT centroids 
with respect to the daughter ground state energies. This 
is again demonstrated in Fig. ^. Here we have chosen to 
plot the centroid energies as function of ( A^ — Z) which is 
the dominating dependence in the parent-daughter mass 
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splitting and is stronger than the Zy{20—Ny) dependence 
of Sp+. Like for GT_, the centroids now group according 
to the pairing structure of the parent ground state. Again 
the centroids are lowest for even-even parents and the 
centroids for the odd- A and odd-odd parents are shifted 
up in energy by about 3 and 6 MeV, resp. 

We note that the same pairing shifts of the 3^+ 
strength has already been suggested by Hansen in a gen- 
eral discussion of the /3 strength in nuclei . 

The dependence of the GT_|_ distribution on the pro- 
ton and neutron number is nicely visualized in fig. ^ for 
the odd-yl = 61 isotone chain. For comparison we note 
that the independent particle model, as assumed by FFN, 
places the GT_|_ centroids at excitation energies of 4 MeV 
for ^-'^Fe and ^^Ni (the daughter nuclei have an even neu- 
tron number) and at 2 MeV for ^^Co. ^^Cu has 9 va- 
lence protons, thus allowing also for a P3/2 proton being 
changed into a pi/2 neutron. Relatedly, the cen- 
troid is shifted slightly to 2.1 MeV. As has already been 
visible in fig. |^, the shell model shows quite a different 
dependence of the GT_|_ distributions. At first, we ob- 
serve that in all cases the strength is concentrated in 
an energy region of about 3-4 MeV width. Further, the 
centroid of this region decreases with increasing neutron 
excess. However, we stress that this decrease is basically 
due to the {N — Z) dependence of the mass difference be- 
tween parent and daughter nucleus, as the GT+ centroid 
energy increases slightly with {N ~ Z) if measured with 
respect to the parent ground state. The figure clearly 
shows no distinct dependence of the GT+ centroid en- 
ergy on the neutron pair configuration. We find the shell 
model centroids at excitation energies of 2.1 MeV (^^Fe), 
3.7 MeV (^^Co), 4.7 MeV (^^Ni), and 6.7 MeV (*^^Cu). 
Clearly the total GT_|_ strength decreases with increasing 
neutron excess in the isotone chain due to the decreas- 
ing number of valence protons and the increasing Pauli 
blocking of the neutrons. 

Finally we remark that the figure also shows that the 
thermal excitation of the strong backresonance transi- 
tions becomes easier with increasing neutron excess, as 
these transitions move to lower excitations energies. 



V. CONCLUSIONS 

Bethe and collaborators had, more than two 
decades ago, focussed the attention on the importance 
played by the weak interaction in a supernova collapse. 
Subsequently about twenty years ago Fuller, Fowler and 
Newman (FFN) ]^-|l0[| outlined in their seminal work the 
theory to calculate stellar weak interaction rates. After 
this pioneering step the problem had been reduced to 
solve the related nuclear structure physics. While FFN 
understood the relevant physics correctly, due to lack of 
data and computational resources they were forced to es- 
timate the relevant weak interaction rates on nuclei in the 
mass range ^ = 45 — 60 phenomenologically. Over the 
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FIG. 9 

A — 61 isotones as function of excitation energy in the daugh- 
ter nucleus. Note the different scales used for the various 
nuclei. 



years experimental findings about the fragmentation and 
positioning of the GT strength in these nuclei indicated 
the need for refinements of the rates and it became obvi- 
ous |jl^;|l^ that the interacting shell model is the method 
of choice for this endeavour. First steps towards this goal 
have been undertaken using the shell model Monte Carlo 
technique ||l^ , but it became clear that shell model diag- 
onalization is the better suited tool to calculate reliable 
stellar rates. Impressive progress in both shell model pro- 
gramming (in particular due to the work by E. Caurier) 
and hardware development allows now for virtually con- 
verged calculations of the GT strength in pf shell nuclei, 
as they play a fundamental role during the presupernova 
collapse. In fact, Caurier et al. have demonstrated |2^] 
that shell model diagonalization is able to reproduce all 
measured GT-(- and GT_ strength distributions on nuclei 
around A « 60 and simultaneously describe the spectra 
and lifetimes of these nuclei also sufficiently well. As this 
is the relevant input to reliably calculate stellar weak 
interaction rates, it has been concluded that the shell 
model diagonalization approach of Ref. has the nec- 
essary predictative power to calculate the rates. In this 
manuscript we have followed this conclusion and have 
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derived stellar weak interaction rates based on state-of- 
the-art shell model diagonalization. The calculation has 
been performed for more than 100 nuclei in the mass 
range A = 45 — 65 and covers the temperature and den- 
sity regime expected in supernova physics. The rates 
have been compiled in a file using the same format as is 
customary for the FFN rates. The electronic version of 
the file can be received from the authors upon request. 
The files are also available in the effective rates formalism 
as derived in [pT| . 

Differences between the shell model and the FFN rates 
are usually related to differences in the placement of the 
GT resonance energies. In analogy to the wellknown sys- 
tematics of the GT_ centroid energies we have derived a 
similar systematics for the GT+ resonance energies. In 
particular we have shown that the GT+ centroid ener- 
gies, if measured with respect to the parent ground state 
energy, are not dependent on the pairing structure of 
the parent nucleus. If the centroids are, however, mea- 
sured with respect to the daughter ground state energy, 
the wellknown pairing energy contributions to the mass 
splitting between parent and daughter nuclei enters and 
the positions of the GT centroids are pushed up in energy 
for odd-A and odd-odd parents by twice and four-times 
the pairing energy, respectively. This systematic has not 
quite been considered in previous estimates of the weak 
interaction rates IH,M,P8^,H • 



EC 

f 



g/cm^. For these nuclei the shell model rates are smaller 
than the FFN rates by factors 39, 12, 10, 346, and 19, 
respectively, where the comparison is made for those tem- 
peratures and densities where Ref. j2^ lists the respec- 
tive nucleus as most important. Already this comparison 
suggests that stellar electron capture rates are noticeably 
smaller than previously assumed. As speculated in [ p^ , 
due to the smaller electron capture rates the core should 
radiate less energy away by neutrino emission, keeping 
the core on a trajectory with higher temperature and 
entropy. 

Electron capture has to compete with (3 decay in the 
stellar environment. For even-even and odd-A nuclei, the 
shell model /3-decay rates are generally also noticeably 
smaller than the FFN rates. But for odd-odd nuclei shell 
model and FFN rates are about the same. As the /3 decay 
of odd-odd nuclei is expected to contribute significantly 
to the total f3 decay rate during a collapse, the total /3- 
decay rate should change less in supernova simulations 
than the electron capture rate, if the FFN compilation is 
replaced by the shell model results. 

Obviously collapse calculations which use the shell 
model stellar weak interaction rates are desired and those 
studies are already initiated. Definite conclusions have 
to wait for the results of these simulations, but we can 
here update one interesting finding put forward in p^ , ^ . 
These authors argued that during the collapse the (3 de- 
cay rate will exceed the electron capture rate for a certain 
range of electron-to-baryon ratios Ye- As a consequence 
the core can radiate energy away, (as the neutrinos can 
still leave the star without interaction) without lowering 
the Ye value. This should have some interesting conse- 
quences for the size of the homologous core. 

In the regime in which the /3-decay rate exceeds 
the electron capture rate has been estimated on the basis 
of shell model rates for a limited set of nuclei. Here we 
update this comparison by adopting the full set of shell 
model rates. To do so, we follow Ref. |44| and define the 
change of due to (3 decay (this increases the charge by 
one unit) and electron capture (which reduces the charge 
by one unit) as 
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FIG. 10. Change in the total electron capture and /3-decay 
rates as a function of electron-to-baryon ratio Ye during a 
supernova collapse. The respective values for the tempera- 
ture and density were taken from the stellar trajectory given 
in H. 



The shell model rates are usually smaller than the FFN 
rates. This is particularly the case for the electron cap- 
ture rate on odd-odd nuclei and to a lesser extent on 
odd-A nuclei. This can be examplified for the nuclei 
55Co,^"Co, ^'^Mn, 6°Co and ^^Mn which the compilation 
of [22] subsequently ranks as the most effective electron 
capturing nuclei in the density regime p ~ 10^ — 10^*^ 



yec(/3) ^ 



dY, 



ec{l3) 



dt 



ec(/3) 
fc 



(20) 



where the sum runs over all nuclear species present in 
the core, and and are the electron capture and /3 
decay rates of nucleus k. The mass fraction is given by 
nuclear statistical equilibrium |p8| . 

As in p3] we will follow the stellar trajectory as given 
in Ref. p^ , although this is expected to change some- 
what if the FFN rates are replaced by the shell model 
rates in the collapse simulation. Fig. |l^ compares 
Y^'^ and Y"/ along the stellar trajectory where Yg re- 
duces here with time. Confirming the results of |Q the 
full set of shell model rates also reveals that the (3 de- 
cay rates are larger than the electron capture rates for 
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Ye = 0.42 — 0.455. This might have important conse- 
quences for the core collapse possibly leading to cooler 
cores and larger Ye values at the formation of the homo- 
loguous core. 

Thielemann and collaborators p5p| have reported first 
attempts to explore the role of the shell model electron 
capture rates in type la supernovae. They find that the 
composition of the matter in the center is less neutron- 
rich than previously assumed. If the FFN rates are re- 
placed by the shell model ones the overproduction of the 
neutron-rich Cr, Ti, and Fe isotopes which has been en- 
countered in previous type la simulations with otherwise 
the same physics input is removed. As the present 
shell model rates likely reduce the uncertainties related 
to the stellar electron capture rates, it is expected that 
type la simulations with the new shell model electron 
capture rates will serve as strict tests for the models and 
their parametrizations. Such comprehensive studies are 
in progress. 
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